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Abstract. The star EZ Peg, long ago classified as cat- 
aclysmic variable, has been shown to be a chromospher- 
ically active binary system of the RS CVn-type. In this 
paper we have analysed, using the spectral subtraction 
technique, simultaneous spectroscopic observations of the 
Ha, H/3, Na i Di and D2, He i D3, Mg i b triplet, Ca 11 
H & K, and Ca 11 infrared triplet lines. We have found 
that the hot component is the active star of the system, 
showing strong emission in the Ha, Ca 11 H & K, He, and 
Ca II IRT lines, and a strong filling-in of the H/3 line, how- 
ever the Na i Di and D2 and Mg i b triplet lines do not 
present filled-in. The He i D3 could present a total filling- 
in due to microflaring. The observed variations (in differ- 
ent epochs and with the orbital phase) of the different ac- 
tivity indicators, formed at different height in the chromo- 
sphere, are correlated. Very broad wings have been found 
in the subtracted profiles of Ha and Ca 11 IRT A8498 and 
A8662 lines. These profiles are well matched using a two- 
component Gaussian fit (narrow and broad) and the broad 
component could be interpreted as arising from microflar- 
ing. The higher luminosity class of the hot component, 
that our spectra seem to indicate, could explain why the 
hot component is the active star of the system. 

Key words: stars: activity - stars: binaries: close - stars: 
chromospheres - stars: late-type - stars: individual: EZ 
Peg 
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the island of La Palma by the Royal Greenwich Observatory 
at the Spanish Observatorio del Roque de Los Muchachos of 
the Instituto de Astrofi'sica de Canarias 



1. Introduction 

EZ Peg (BD -1-24° 4742) has a misleading history of con- 
fused interpretation of the origin and composition of the 
system. The spectrum of EZ Peg was originally classified 
by Cannon as G5 (Schlesinger et al. 1934). Vyssotsky & 
Balz (see Alden 1958) claimed that the star varies from BO 
to G5, due to a spectrum taken in 1943 which appeared 
to be a B star spectrum. Later Howell & Bopp (1985) pre- 
sented evidence that this erroneous interpretation of the 
spectrum was due to poor resolution data. Irvine (1972) 
proposed that EZ Peg was an U Gem system (a cata- 
clysmic variable), but it is an unlikely explanation due to 
the normal UV emission of the system: in a color-color 
plot there is no strong UV excess typical of the dwarf 
novae and symbiotics (Szkody 1977). Furthermore linear 
polarimetry measurements of the system (Szkody et al. 
1982) show that EZ Peg does not seem to be a cataclysmic 
variable. 

Howell & Bopp (1985) observed the double-hned spec- 
trum of EZ Peg and interpreted the color-indices (B-V) 
and (U-B) as that of a composite spectral type of G5 IV 
+ KG IV. They observed Ca 11 H & K, Ha and He hues 
in emission, and they supposed that this emission origi- 
nated in the cooler star, because it is the most typical case. 
Howell et al. (1986) suggested that both components of EZ 
Peg were subgiants, and pointed out a slightly asymmet- 
ric proflle in Ha, interpreted as a vertical motion in the 
atmosphere of the active star of the system. 

Griffin (1985) presented radial- velocity observations of 
the system, and he calculated the orbital elements. Re- 
cently Gunn et al. (1996) made cross-correlation radial ve- 
locity measurements of EZ Peg, showing a disagreement 
with the phases calculated from the ephemeris given by 
Griffin (1985). Vilkki et al. (1986) calculated a distance to 
EZ Peg of 83 pc using the trigonometric parallax method, 
and they suspected that EZ Peg belongs to the Hyades 
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cluster. Barret (1996) calculated a distance of 82 pc using 
linear polarimetry. 

Finally, recent studies in the IR band by Mitrou et al. 
(1996) with the IRAS satellite indicated that the system 
does not seem to present any circumstellar matter. 

Continuing with the series of multiwavelength optical 
observations of chromospherically active binaries (Montes 
et al. 1997b, hereafter Paper I), in this paper we present 
simultaneous observations in different optical wavelength 
bands, including also a high resolution echelle spectrum. 
All the data seem to point out that EZ Peg is a RS CVn 
system with an evolved G5 active star, and a KOV or IV 
non-active star. We have applied the spectral subtraction 
technique, observing the Ca ii H & K, He, Ha and Ca IR 
triplet lines in emission, and a filling-in of the H/? line, 
proceeding from the hotter component of the system. 

In Sect. 2 we give the details of our observations and 
data reduction. In Sect. 3 we discuss some of the stellar 
parameters of the binary system and in Sect. 4 the be- 
havior of the different chromospheric activity indicators is 
described. Finally in Sect. 5 we present the conclusions. 

2. Observations and Data Reduction 

Spectroscopic observations of EZ Peg in several optical 
chromospheric activity indicators have been obtained dur- 
ing two observing runs with the 2.5 m Isaac Newton Tele- 
scope (INT) at the Observatorio del Roque de Los Mucha- 
chos (La Palma, Spain) using the Intermediate Disper- 
sion Spectrograph (IDS). In the first run (September IS- 
IS 1995) we have obtained two spectra: in the Ha and 
Na I Di, D2, He I D3 line regions, using grating H1800V, 
camera 500 and a 1024 x 1024 pixel TEK3 CCD as de- 
tector. The reciprocal dispersion achieved is 0.24 A/pixel 
which yields a spectral resolution of 0.48 A and a useful 
wavelength range of 250 A centered at 6563 A (Ha) and 
5876 A (Hei D3) respectively. In the second run (Novem- 
ber 26-28 1996) we used the same configuration but two 
different gratings H1200R and H1200B, for the red and 
blue spectral regions respectively. Using the same config- 
uration a reciprocal dispersion of 0.39 A was obtained 
which yields a spectral resolution of 0.94 A and a use- 
ful wavelength range of 385 A centered at 6563 A (Ha), 
5876 A (Hei D3), 4861 A (H/3), 3950A (Ca 11 H&K) 
respectively. 

In Table |l| we give the observing log. For each spectral 
region observed we hst the date, UT, orbital phase {^p) 
and signal to noise ratio (S/N) obtained. 

We also analyse here echelle spectra obtained with 
the 4.2 m WiUiam Herschel Telescope (WHT) and the 
Utrecht Echelle Spectrograph (UES) on July 31 1993 that 
we have retrieved from La Palma Data Archive (Zuider- 
wijk et al. 1994). These WHT/UES spectra were obtained 
with echelle 31 (31.6 grooves per mm) and a 1280 x 1180 
pixel EEV6 CCD as detector. The central wavelength 
is 6127 A and covers a wavelength range from 4845 to 



8805 A over 53 echelle orders. The achieved reciprocal 
dispersion ranges from 0.045 to 0.080 A/pixel. In order to 
improve the S/N ratio and to avoid cosmic ray events, the 
final spectrum analysed here is the resulting from com- 
bine 28 CCD images with exposure times of 240 s each, 
obtaining a final spectrum with a S/N of 315 in the Ha 
line region. In Table § we give the wavelength range and 
the spectral lines of interest in each echelle order. 

In addition we also use one spectrum in the region 
of the Ca 11 IRT lines taken on January 16 1997 with 
the 1.0 m Jacobus Kapteyn Telescope (JKT) using the 
Richardson-Brealey Spectrograph (RBS) with a grating 
of 1200 grooves per mm, and a 1024 x 1024 pixel TEK4 
CCD as detector. The reciprocal dispersion achieved is 
0.85 A/pixel which yields a spectral resolution of 1.7 A 
and a useful wavelength range of 865 A centered at 
8765 A. 

The spectra have been extracted using the stan- 
dard reduction procedures in the IRAF El pack- 
age (bias subtraction, flat-field division, and opti- 
mal extraction of the spectra). The reduction of the 
JKT/RBS spectrum was carried out using the reduc- 
tion package ReD^E (Cardiel & Gorgas 1997, see also 
http: //www. ucm.cs/OTROS/Astrof/reduceme/reduceme. html). 
The wavelength calibration was obtained by taking spec- 
tra of a Cu-Ar lamp in the INT/IDS spectra, a Th-Ar 
lamp in the WHT/UES spectrum and a Cu-Ne lamp in 
the JKT/RBS spectrum. Finally, the spectra have been 
normalized by a polynomial fit to the observed continuum. 

3. Stellar parameters of the binary system EZ Peg 

EZ Peg is a double-lined spectroscopic binary with an or- 
bital period of 11.6598 days (Grifhn 1985), classified as 
G5V-IV/K0IV. It is in synchronous rotation with a pho- 
tometric period of 11.6626 days (Howell et al. 1986). In Ta- 
ble H we show the adopted stellar parameters from Strass- 
meier et al. (1993) except for Tconj (see below). 

3.1. Radial velocities 

In our high resolution echelle spectrum we can resolve 
both components in all the orders. We have calculated 
the radial velocity difference for this observation using 39 
lines clearly seen in the echelle spectrum as double lines. 
The value obtained is 16.17±0.72 km s~^ which results in 
a value for the orbital phase of 0.944±0.002. 

According to the ephemeris given by Strassmeier et al. 
(1993) the calculated phase for this spectrum is 0.997, 
which is clearly in disagreement with our observation, 
since at this phase the radial velocity difference is 0.92 km 
s~^ which yields a difference in wavelength between both 

^ IRAF is distributed by the National Optical Observatory, 
which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under contract with the National 
Science Foundation. 
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Date IDS-INT UES-WHT 



Na I Di, D2, He I D3 Ca II H fc K H/3 Echcllc 



UT 1^ S/N UT 1/7 S/N UT S/N UT ip S/N UT S/N 

1993/07/31 - _ - _ _ _ _ . - - - _ 02:05 0.954 315 

1995/09/14 00:36 0.416 295 00:04 0.414 362 _ . . ... _ _ 

1995/09/15 00:23 0.501 332 00:02 0.500 343 _ _ _ _ _ _ _ _ 

1996/11/26 21:25 0.141 295 _ _ _ 23:11 0.148 295 _ _ _ _ _ 

1996/11/27 20:39 0.224 295 20:49 0.225 362 19:35 0.220 295 _ _ _ . . 

1996/11/28 19:24 0.306 295 20:25 0.310 362 21:44 0.314 295 21:19 0.312 362 



Table 3. Stellar parameters 



Tap SB R d B-V V-R Tt-onj P^rb Prot Vsini 

(Ro) (pc) (H.J.D.) (days) (days) (kms-^) 

G5V-IV/K0IV: 2 - 83 0.73/0.91 0.58/0.71 2445737.165 11.6598 11.6626 9/7 



components of only 0.02 A. This AA is lower than the reso- 
lution of the spectrum and then the lines from both com- 
ponents should be blended. Similar problems have been 
found in the other spectra when we use this ephemeris. 
The phase deviation is the same for all the observations 
and is very large to be ascribed to orbital period varia- 
tions, so we thought of an error in the heliocentric Julian 
date on conjunction (Tconj)- In effect, we have found a 
mistake in the Tconj given in the Strassmeier et al. (1993) 
catalog where the orbital data by Griffin (1985) are com- 
piled. In the original paper the epoch is determined using 
Modified Julian Date (MJD) and so a difference of 0.5 
days must be added to the Strassmeier et al. (1993) date 
which used 2440000.0 Juhan day as reference date. The 
orbital phase of the echelle spectrum obtained using the 
corrected value of Tconj (2445737.165) is 0.954. Now, the 
difference with the orbital phase that we deduce from our 
spectrum is only 0.01, which is within the uncertainties 
given by Griffin (1985). We would notice that the large 
discrepancies in the radial velocities which Gunn et al. 
(1996) found could be reduced using the corrected Tconj. 
Therefore we have used the corrected Tconj given in Ta- 
ble H to calculate the orbital phases of all the observations 
reported in this paper. 

3.2. Spectral type and luminosity class 

Several Ti I lines can be used as a luminosity classification 
criterium (Keenan & Hynek 1945; Kirkpatrick et al. 1991; 
Ginestet et al. 1994; Jashek & Jashek 1995) because they 
present a positive luminosity effect (i.e. they are stronger 
in giants than in dwarfs). We have analysed the behaviour 
of the Ti I lines that appear in our echelle spectrum in 



order to improve the luminosity classification of the EZ 
Peg components. 

In Fig. ^ we can notice the enhancement of the multi- 
plet 33 Ti I lines for the system in comparison with main 
sequence stars of spectral type similar to the components 
of EZ Peg. It is worth mentioning the case of the Ti I 
A8382 line which is even more intense than the Fe I A8388 
line, showing a different behaviour that in the subgiant 
reference star (see Fig. ^ lower panel). This behaviour 
can be also observed in other Ti I lines like A5366.651 
(Multiplet 35), A5446.593 (M. 3), A5866.453, A5937.806 
and A5941.755 (M. 72), A6554.226 and A6556.066 (M.102), 
A6599.112 (M.19). These results suggest a clear difference 
in luminosity class between the two components of the 
system, being the G5 primary star more evolved than a 
subgiant, meanwhile the KO star is probably a luminosity 
class V star. 

As is explained in the following section all the activ- 
ity indicators analysed here show that the hot compo- 
nent is the active star of the system, contrary to the usual 
behaviour observed in chromospherically active binaries. 
Only in some BY Dra stars the hot component tends to 
be the most active star of the system or even the only 
active component (Montes et al. 1996a). So we have used 
our Ca II K spectra to determine the absolute visual mag- 
nitude, Mv of the active component by the application 
of the Wilson-Bappu eflFect (Wilson & Bappu 1957). The 
mean emission line width measured in our Ca ii K spec- 
tra is 1.229 A, which come out to 60.39 km s^^ after the 
quadratic correction of the instrumental profile. We have 
used the relation found for chromospherically active bina- 
ries by Montes et al. (1994) obtaining My = 2.27. This 
value is lower than the My that corresponds to a G5V 




Fig. 1. WHT/UES (Jul 1993), spectra in the region of several lines of Ti I (M.33). In the upper panel we plot EZ Peg and in 
the lower panel three reference stars: HD 186427 (G3V), HD 44867 (G8IV), HD 185144 (KOV). 
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Table 2. UES spectral orders 



No. 


X- 




Chromospheric Lines 


Other lines 


1 


4845.80 


4897.29 


H/3 






2 


4887.53 


4939.46 








3 


4929.98 


4982.35 








4 


4973.16 


5026.00 








5 


5017.11 


5070.41 








6 


5061.85 


5115.62 








7 


5107.38 


5161.63 








8 


5153.74 


5208.48 


Mg i b A5167-72-83 






9 


5200.94 


5256.19 








10 


5249.02 


5304.77 








11 


5297.99 


5354.27 








12 


5347.89 


5404.69 




Ti i 


L A5366 


13 


5398.73 


5456.07 




Ti i 


L A5426 


14 


5450.54 


5508.43 








15 


5503.36 


5561.81 








16 


5557.21 


5616.23 








17 


5612.13 


5671.73 








18 


5668.14 


5728.34 








19 


5725.28 


5786.09 








20 


5783.58 


5845.01 








21 


5843.08 


5905.14 


Na i Di, D2, He i D3 


Ti i 


L A5866 


22 


5903.82 


5966.53 








23 


5965.84 


6029.20 








24 


6029.17 


6093.21 








25 


6093.87 


6158.59 








26 


6159.96 


6225.39 








27 


6227.51 


6293.66 








28 


6296.55 


6363.44 








29 


6367.15 


6434.78 




Fe i 


L A6411.66 


30 


6439.35 


6507.75 




Fe i 


L A6495 


31 


6513.20 


6582.39 


Hq 






32 


6588.77 


6658.76 








33 


6666.11 


6736.93 




Fe i 


L A6678 










Li i 


A6708 










Ca 


i A6718 


34 


6745.30 


6816.95 








35 


6826.39 


6899.90 




O2 


A6867 


36 


6909.45 


6982.85 








37 


6994.56 


7068.87 








38 


7081.80 


7157.03 








39 


7171.24 


7247.43 








40 


7262.98 


7340.14 








41 


7357.09 


7435.25 








42 


7453.68 


7532.86 








43 


7552.84 


7633.08 




B-band A7600 


44 


7654.68 


7736.00 




K i 


A7665, 7699 


45 


7759.30 


7841.73 








46 


7866.82 


7950.40 








47 


7977.38 


8062.12 








48 


8091.08 


8177.03 








49 


8208.08 


8295.27 








50 


8328.51 


8416.99 




Ti i 


L (M.33) 


51 


8452.54 


8542.33 


Ca ii IRT A8498, 8542 


Ti i 


L A8468 










Fe : 


L A8514 










Ti i 


L A8518 


52 


8580.32 


8671.46 


Ca ii IRT A8662 


Fe ! 


L A8621 


53 


8712.02 


8804.56 









(Mv(Tsp) = 5.1 from Landolt-Bornstcin (Schmidt-Kaler 
1982)) but higher than for a G5III (Mv(Tsp) = 0.9). This 
result indicates that the active component of EZ Peg may 
be of luminosity class IV or higher. 

The difference in luminosity class between the hot and 
cool components that our spectra indicate could explain 
why the hot component is the active star (see below) , since 
if the hot component is a more evolved star it could have 
developed a deep convective zone and therefore, according 



to the dynamo mechanism, should present a higher activ- 
ity level. Anyway the cooler KOV component could also be 
active and present a faint emission in the H & K lines, but 
due to the small contribution of that star to the combined 
spectrum of the system, this emission is very difficult to 
be observed. 

3.3. The Lii X6707.8 line 

We also analyse here the Lithium abundance of EZ Peg 
using the Li i A6707.8 line. This line appears in the order 
33 of the WHT/UES echelle spectrum together with the 
Fe I A6678 and Ca i A6718 fines. In Fig. || we can clearly 
see the Fe i and Ca i lines from both components, the 
expected positions of these features for the hot (H) and 
cool (C) are given in this figure. However, the absorp- 
tion Li I line appears centered in the wavelength-position 
corresponding to the hot component, so we can suppose 
that the absorption line comes only from this star or that 
the contribution of the cool component to the observed 
Li I is very small. This assumption seems to be reason- 
able because, as we show below, the hot component is 
the active star of the system and a high Li i abundance 
is found in chromospherically active binaries (Fernandez- 
Figueroa et al. 1993; Barrado et al. 1997). In order to ob- 
tain the Li I equivalent width (EW), we have corrected the 
total EW measured, EW(Li i-fFe i)=50.4 (niA) from the 
relative contribution to the continuum (0.65). Then the 
contribution from the blend with the Fe i A 6704.41 line 
must be subtracted. The EW of the Fe i line was calcu- 
lated from the empirical relationship with (B-V) given by 
Soderblom et al. (1990), obtaining EW(Fe i) = 15.8 (mA), 
the same value is found with the relationship given by Fa- 
vata et al. (1993) using the effective temperature. Finally, 
the corrected EW(Li i)corr=61.7 (mA) was converted into 
abundances by means of the curves of growth computed 
by Pallavicini et al. (1987), obtaining log N(Li i)=2.15 
(on a scale where log N(H)=:12.0) with an accuracy of 
w 0.30 dex. This high Li i abundance is similar to that 
obtained for other chromospherically active binaries by 
Fernandez-Figueroa et al. (1993) and Barrado et al. (1997) 
and it is related to the high activity level of EZ Peg. 

4. Chromospheric activity indicators 

In the following we describe the behaviour of the differ- 
ent optical chromospheric activity indicators (formed at 
different atmospheric heights) observed for EZ Peg: Na I 
Di, D2, and Mg l b triplet (upper photosphere and lower 
chromosphere), Ca 11 IRT lines (lower chromosphere). Ha, 
H/3, Can H & K (middle chromosphere), and He i D3 
(upper chromosphere). The chromospheric contribution in 
these features has been determined using the spectral sub- 
traction technique described in detail by Montes et al. 
(1995a, b, c), (see also Paper I). The synthesized spec- 
trum has been constructed with reference stars of spectral 
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Fig. 2. WHT/UES (Jul 1993), spectrum in the region of the Li I 6707.8 A line including also the Fe I 6677.99 A and Ca I 
6717.69 A . The expected positions of this features for the hot (H) and cool (C) components are given. 



types G5V and KOIV, taken from the spectral library of 
Montes et al. (1997a), in the case of the INT/IDS spec- 
tra. For the WHT/UES spectrum we have used a G3V 
and KOV taken during the same observing run, and for 
the JKT/RBS spectrum a G5III-IV and KOV have been 
used. The use of reference spectra of spectral type slightly 
different for each observing run and different also from 
the system component spectra introduce some residual in 
the absorption lines, but it does not affect significatively 
the excess emission equivalent width measures in the sub- 
tracted spectra. The uncertainties are within similar or- 
der to that obtained in the reduction, synthesis and fitting 
procedures. Furthermore, we have considered a clear de- 
tection of excess emission or absorption lines only when 
these features in the difference spectrum are larger than 3 
cr. 

The contribution of the hot and cool components to 
the total continuum that yields a better fit between the 
observed and subtracted spectra is 0.65/0.35. 

The line profiles are displayed in Figs. 1-8. For each 
spectrum we plot the observed (solid-line) , the synthesized 
(dashed-hne) and the subtracted one, additively offset for 
better display (dotted hne). 

4-.1. The Ha line 

We have taken several spectra of EZ Peg in the Ha line 
region in three different epochs. One WHT/UES spec- 
trum in July 1993 at the orbital phase 0.954 (Fig. |), 
two INT/IDS spectra in September 1995 at the phases 
0.416 and 0.501 (Fig. |) and three more INT/IDS spectra 
in November 1996 at the phases 0.141, 0.224 and 0.306 
(Fig. ^). In all the spectra we can see the Ha line in emis- 
sion above the continuum with an asymmetry in the wave- 
length position corresponding to the Ha absorption line of 
the cool component. The spectral subtraction reveals that 
the hot star is responsible for the excess Ha emission. 



Table I gives the Ha line parameters, measured in the 
observed and subtracted spectra of EZ Peg as follows: col- 
umn (1) the observing run, column (2) the orbital phase 
{ip) and columns (3), (4), (5) the following parameters 
measured in the observed spectrum: the full width at half 
maximum (Wobs); the residual intensity, Rc; and the Ha 
core flux, F(1.7A), measured as the residual area below 
the central 1.7 A passband. The last four columns give 
the following parameters measured in the subtracted spec- 
trum: the full width at half maximum (Wgub), the peak 
emission intensity (I), the excess Ha emission equivalent 
width (EW( Ha)), and absolute fluxes at the stellar sur- 
face logFg (Ha) obtained with the calibration of Pasquini 
& Pallavicini (1991) as a function of (V - R), very simi- 
lar values of Eg (Ha) are obtained using the more recent 
calibration of Hall (1996) as a function of (V - R) and (B - 
V). As can be seen in Table ^ and Figures ^, ^ and ^ the 
excess Ha emission of EZ Peg show small variations with 
the orbital phase and also seasonal variations from Jul 93 
to Nov 96. The higher Ha EW has been reached in the 
first night of the 1996 run and in the 1993 observation. 
The higher values correspond to the orbital phases closer 
to the conjunction (0.0). 

The Ha subtracted profiles present at all the epochs 
broad wings which are not well matched using a single- 
Gaussian fit. These profiles have therefore been fitted 
using two Gaussian components. The parameters of the 
broad and narrow components used in the two-Gaussian 
fit are given in table ^ and the corresponding profiles are 
plotted in the left panel of Figures 1 1 and |. In some cases 
the blue wing is noticeable stronger than the red wing and 
the fit is better matched when the broad component is 
blue-shifted with respect to the narrow component. The 
cool component of EZ Peg cannot account for this effect 
since any emission from that star would be on the red side 
of the observed line for these orbital phases. We have also 
observed broad components in the most active systems of 
the sample of Paper I which we interpreted as microflar- 
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Observed Ha Spectrum Subtracted Ha Spectrum 



Obs. 


V 


Wobs 


Rc 


F(1.7A) 




I 


EW 


log Fg 






(A) 






(A) 




(A) 




UES 1993 


0.954 


0.96 


1.325 


1.987 


1.62 


0.882 


2.085 


7.06 


INT 1995 


0.416 


0.65 


1.139 


1.695 


1.70 


0.655 


1.488 


6.91 




0.501 


0.54 


1.059 


1.657 


1.69 


0.639 


1.400 


6.89 


INT 1996 


0.141 


1.48 


1.131 


1.779 


2.43 


0.604 


2.272 


7.10 




0.224 


1.55 


1.157 


1.812 


2.30 


0.634 


1.862 


7.01 




0.306 


1.35 


1.134 


1.784 


2.16 


0.638 


1.677 


6.96 



Table 5. Parameters of the broad and narrow Gaussian components used in the fit of the Hq subtracted spectra 









Ha broad component 
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UES 1993 


0.954 


0.273 


3.935 


1.143 


54.8 


0.648 


1.367 


0.942 


45.2 


INT 1995 


0.416 


0.152 


4.035 


0.651 


43.7 


0.535 


1.470 


0.837 


56.3 




0.501 


0.129 


4.096 


0.562 


40.1 


0.535 


1.470 


0.838 


59.9 


INT 1996 


0.141 


0.123 


9.418 


1.209 


53.2 


0.493 


2.028 


1.064 


46.8 




0.224 


0.124 


5.622 


0.740 


39.7 


0.520 
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1.122 


60.3 




0.306 
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0.536 


32.0 


0.572 


1.876 
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Fig. 3. WHT/UES Ha spectrum (Jul 1993). In the left panel we plot the observed spectrum (solid-line), the synthesized 
spectrum (dashed-line), the subtracted spectrum, additively offset for better display (dotted line). In the right panel we plot 
the subtracted Ha profile (dotted line), and the the two Gaussian components fit (solid-line). The sort-dashed-line represents 
the broad component and the large-dashed-line the narrow one. 



ing that occurs in the chromosphere by similarity with the 
broad components also found in the chromospheric Mg ii 
h & k lines (Wood et al. 1996) and in several transition 
region lines of active stars (Linsky & Wood 1994; Linsky 
et al. 1995; Wood et al. 1996, 1997; Dempsey et al. 1996a, 
b; Robinson et al. 1996). The contribution of the broad 
component to the total EW of the line of EZ Peg ranges 
between 32% and 55% which is in the range observed to 
the stars analysed in Paper I. The larger changes in the 



excess Ha emission appear to occur predominantly in the 
broad component and its contribution is related to the de- 
gree of stellar activity. It is noticeable the change in the 
wings of the Ha line from the 1st to the 2nd night of the 
1996 run, the two Gaussian components fit reveals that 
the contribution of the broad component changes from 
53.2% to 39.7% (see Fig|). 
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Fig. 4. INT/IDS Ha spectra (Sep 1995), as in Fig || 



4.2. The HP line 

Two spectra in the H/3 line region are available. The 
WHT/UES spectrum (July 1993) and other in the 
INT/IDS 1996 run. In both spectra the H/3 line appears 
in absorption and the application of the spectral subtrac- 
tion technique reveals a clear excess emission in this line 
(see Fig. ^ . Again the relative velocity shift shows that 
the emissions belong to the hotter component. Contrary 
to the case of the Ha line the H/3 subtracted profiles do 
not present broad wings and they are well matched using 
a single- Gaussian fit. 

Table ^ gives the H/3 line parameters, measured in the 
observed and subtracted spectra as in the case of the Ha 
line, see also our previous study of the excess H/3 emission 
in active binaries (Montes et al. 1995d). 

In Table |6| we also give the ratio of excess emission 



EW, 



EW(Hq) 



for the two H/3 observations of EZ Peg, and 



EW(H/3) ' 

the ratio of excess emission with the correction: 



E 



Ho 



0.2444 * 2.512(^-«). 
Eh/3 EW(H/3) 

given by Hall & Ramsey (1992) that takes into account the 
absolute flux density in these lines and the color difference 
in the components. 

We have used this ratio as a diagnostic for discrimi- 
nating between the presence of plages and prominences 
in the stellar surface, following the results of Hall & 



Ramsey (1992) who found that low EHa/En/s (« 1-2) 
can be achieved both in plages and prominences viewed 
against the disk, but that high ratios (k, 3-15) can only 
be achieved in extended regions viewed off the limb. The 
high ratio (EHa/En/s >3, see Table ^ that we have found 
in EZ Peg indicates that the emission would arise from 
extended regions (prominences). 

The excess H/3 emission is somewhat larger in the 1993 
observation than in the last night observation of the 1996 
run, in agreement with the behaviour observed for the Ha 
line. 



4.3. The Na I Di and D2 and He I Dt, lines 

We have five spectra of EZ Peg in the region of Na i Di 
and D2 and He i D3 lines: one WHT/UES spectrum in 
July 1993 at the orbital phase 0.954 (Fig. |), two INT/IDS 
spectra in September 1995 at the phases 0.414 and 0.500 
and two spectra in INT/IDS November 1996 at the phases 
0.225 and 0.310 (Fig. |). We have not detected any signif- 
icant absorption or emission for He i D3. 

Saar et al. (1997) in their study of He i D3 in G and 
K dwarfs found that while for few active stars the He I 
D3 line behaves " normally" , increasing in absorption with 
increasing rotation, and showing consistent correlations 
with other activity indicators, the behavior clearly di- 
verges when stars become very active. They found a large 
filling-in He i D3 line in some very active K stars. On 
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Fig. 5. INT/IDS Ha spectra (Nov 1996), as in Fig 
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Fig. 6. WHT/UES H/3 (Jul 1993) and INT/IDS H/? (Nov 1996) spectra. For each spectrum we plot the observed spectrum 
(solid- line) , the synthesized spectrum (dashed-line), the subtracted spectrum, additively offset for better display (dotted line), 
and the Gaussian fit to the subtracted spectrum (dotted-dashed line) 
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Table 6. H/3 line measurements in the observed and subtracted spectrum 



Observed H/3 Speetrum Subtracted H/3 Spectrum 

Obs. V W„bs Rc F(1.7A) W,,b I EW 

(A) (A) (A) 

UES 1993 0.954 0.91 0.561 1.235 0.94 0.364 0.489 4.27 3.49 

INT 1996 0.312 2.13 0.660 1.237 1.54 0.244 0.411 4.08 3.33 



the other hand, they also found He i D3 EW significantly 
lower than the theoretical maximum value given by An- 
dretta & Giampapa (1995), suggesting that the line could 
be filling-in due to frequent low-level flaring. 

Both facts could explain the absence of emission or ab- 
sorption due to a total filling-in of the He I D3 line that 
we have found in EZ Peg and also in other dwarfs and 
subgiants chromospherically active binaries (see Paper I), 
since these stars are very active and also present low- 
level flaring (microflares) as indicated by the Ha broad 
component that we have found (see sect. 4.1 and Pa- 
per I). In some cases He i D3 goes into emission due to 
flares (Huenemoerder & Ramsey 1987; Huenemoerder et 
al. 1990; Montes et al. 1996b; Paper I). However, in the 
most evolved stars the behaviour is different as conse- 
quence of the lower chromospheric densities of these stars, 
as it is suggested by Saar et al. (1997) and as it is indi- 
cated by the presence of the He i D3 in absorption that 
we have reported in Paper I for several chromospherically 
active binaries with giant components. 

In the Na i Di and D2 lines the spectra are too noisy 
to distinguish a clear flUed-in. Once again we can observe 
(see Fig. ^ that the redward lines of the double- lined spec- 
trum usually are the deepest, but in the Na i D whose 
blueshifted components are the deepest ones. As it is well 
known, the Na i D lines are very sensitive to the tempera- 
ture; this behaviour of the lines confirms that the hot star 
of the system is the redshifted one, which, as mentioned 
above, shows very high activity in Ha. 

4-. 4- The Mg i h triplet lines 

The Mg I b triplet lines A5167.33, A5172.68, and A5183.61 
are formed in the lower chromosphere and the region of 
temperature minimum and they are good diagnostics of 
photospheric activity (Basri et al. 1989). A small flUing-in 
of these lines have been recently found in some RS CVn 
system (Gunn & Doyle 1997; Gunn et al. 1997). As in the 
case of the Na i Di and D2 lines, in the subtracted spec- 
trum corresponding to the echelle images no clear fiUed-in 
has been found. Thus the photospheric activity in EZ Peg 
is very low. 



4.5. The Call H &z K and He lines 

We have taken three spectra in the Ca 11 H & K line region 
during the 1996 run (see Fig. These spectra exhibit 
strong Ca 11 H & K emission lines and a small bump in 
the wavelength-position corresponding to the He line. The 
H & K emission lines appear blue-shifted with respect to 
the observed absorption lines in agreement with the or- 
bital phase, which conflrms that the emission originates 
in the hot component. The spectral subtraction yields a 
satisfactory flt and reveals a small but clear emission in 
the He line. 

The excess Ca 11 H & K and He emissions equivalent 
widths (EW) have been measured by reconstruction of the 
absorption line proflle (described by Fernandez-Figueroa 
et al. 1994) and using the spectral subtraction technique 
(explained by Montes et al. 1995c, 1996a). The Ca 11 H 
& K surface flux, Fs(Ca 11 H & K) has been obtained by 
means of the linear relationship between the absolute sur- 
face flux at 3950 A (in erg cm^^ s^^ A^^) and the colour 
index (V-R) by Pasquini et al. (1988). 

Table gives the Ca 11 H & K and He line parameters, 
measured in the observed and subtracted spectra. Column 
(1) gives the observing run. Column (2) gives the orbital 
phase {ip). In columns (3) and (4) we list the EW for the K 
and H lines, obtained by reconstruction of the absorption 
line profile, in columns (5) to (7) we give the EW for the 
K, H and He lines, measured in the subtracted spectrum, 
and in columns (8) to (10) the corresponding surface fluxes 
are listed. 

The excess Ca 11 H & K and He emissions change with 
the orbital phase during the three nights of the 1996 run in 
the same way as the corresponding excess Ha emissions. 

4-6. The Call IRT lines 

The Ca 11 infrared triplet (IRT) A8498, A8542, and A8662 
lines are other important chromospheric activity indica- 
tors (Linsky et al. 1979; Foing et al 1989; Dempsey et al. 
1993). These lines share the upper levels of the Ca 11 H 
& K transitions (see for example Mihalas 1978), but un- 
like the blue Ca 11 lines, the IRT lies in a region of the 
spectrum with a well-defined continuum, making the cal- 
ibration simpler. The IRT is formed in the lower chromo- 
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Fig. 7. WHT/UES (Jul 1993), spectrum in the region of the NaT and Hel hues including also the Cal 5857.454 A, Fel 5862.357 A 
and Til 5866. 453A. The expected positions of these features for the hot (H) and cool (C) components are given. 
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Fig. 8. INT/IDS spectra in the region of the Nal and Hel lines (Sep 1995) and (Nov 1996). For each spectrum we plot the 
observed spectrum (solid- line), the synthesized spectrum (dashed- line) , the subtracted spectrum, additively offset for better 
display (dotted line). 

Table 7. Ca 11 H &: K and He lines measures in the observed and subtracted spectrum 
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6.91 


6.79 6.28 
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0.671 


1.191 
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Fig. 9. INT/IDS Ca II H & K spectra (Nov 1996), as in Fig .| 



sphere. Shine & Linsky (1972, 1974) observed solar IRT 
plage profiles concluding that the emission results from 
a steepening of the temperature gradient in the chromo- 
sphere. This forces the temperature, electron density, and 
Ca II density population to be higher for a given optical 
depth, compared to the quiet chromosphere, thus increas- 
ing the source function thereby producing emission in the 
line cores. 

In the observations carried out by Dempsey et al. 
(1993) only one star of the sample (BY Dra) showed emis- 
sion above the local continuum in only the A8498 line, and 
other seven stars revealed emission rising above the shal- 
low line core. 



The A8498 and A8662 fines are included in the 1993 
WHT/UES echelle spectrum of EZ Peg (see Fig. |o|), and 
both lines exhibit a strong emission superposed to the red 
wing of the corresponding absorption line, in agreement 
with the orbital phase (0.954). In the 1997 JKT/RBS spec- 
trum at orbital phase 0.511 the three Ca ii IRT lines are 
included (see Fig. |ll|). In this case the lines show a small 
emission reversal centered at the absorption line in the 
observed spectrum and a clear excess emission in the dif- 
ference spectrum. 

Table || gives the Ca II IRT line parameters, measured 
in the observed spectra by reconstruction of the absorp- 
tion fine profile, and using the spectral subtraction. The 
columns in this tabic have the same meaning as for the 
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Fig. 10. WHT/UES Ca ll IRT A 8498 and A 8662 spectra (Jul 1993), as in Fig 
Ca II IRT (JKT/RBS Jan 1997) 
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Fig. 11. JKT/RBS Ca ll IRT spectrum (Jan 1997) at orbital phase 0.511 



Ca II H & K lines (see Table 0) . Notice that the under- 
estimation of the reconstruction method is larger for the 
IRT lines than for the Ca ii H & K lines. The absolute 
fluxes at the stellar surface have been obtained using the 
calibration of Hall (1996) as a function of (V - R). 

As in the case of the Ha line the subtracted profiles 
of the A8498 and A8662 lines also present broad wings. 
Therefore we have used, as for the Ha profile, a two gaus- 



sian (broad and narrow) fit. In Table ^ we list the param- 
eters (I, FWHM, EW) of the broad and narrow compo- 
nents. The narrow components of both lines are very sim- 
ilar; however, the contribution of the broad component to 
the total EW is larger in the A8662 line. In the subtracted 
spectrum the intensity of the A8498 line is higher than 
in A8662; however, due to the contribution of the broad 
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Table 8. Ca ii IRT lines measured in the observed and subtracted spectra 
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UES 1993 
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Table 9. Parameters of the broad and narrow Gaussian components used in the fit of the Ca ii IRT subtracted spectra 
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FWHM 
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EWb/EWt 
(%) 
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FWHM 
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Ca ii IRT line A8498 


UES 1993 


0.954 


0.089 


2.497 


0.237 


34.8 


0.693 


0.600 


0.442 


65.2 


Ca ii IRT line A8662 


UES 1993 


0.954 


0.086 


5.393 


0.452 


50.7 


0.564 


0.733 


0.440 


49.3 



component, the total excess emission EW is larger in the 
A8662 line. 

5. Conclusions 

In this paper we have analysed, using the spectral subtrac- 
tion technique, simultaneous spectroscopic observations of 
several optical chromospheric activity indicators, at differ- 
ent epochs, of the RS CVn binary system EZ Peg. 

We have found the Ha line in emission above the con- 
tinuum in all the spectra, with noticeable changes in the 
line profile, that in part are due to the superposition of 
the absorption line to the other component. The sub- 
tracted Ha profile is well matched using a two-components 
Gaussian fit (narrow and broad). The broad component 
is mainly responsible for the observed variations of the 
subtracted profile, and its contribution to the total EW 
increases with the degree of activity. We have interpreted 
this broad component (also observed in other active sys- 
tems in Paper I) as arising from microflaring activity that 
take place in the chromosphere. 

A filled-in in the H/3 line has been found; however, the 
subtracted profiles do not present broad wings. The ratio 
(EHa/En/?) that we have found indicates that the emission 
would arise from extended regions viewed off the limb. 

Strong Ca Ii H & K emission from the hot component 
is observed in our spectra. The spectral subtraction re- 
veals also a small but clear emission in the He line. The 
fines A8498 and A8662 of the Ca il IRT also show a strong 
emission. As in the case of the Ha line the subtracted pro- 
files of these lines are well matched using a two-component 
Gaussian fit. 



We suggest that the He I D3 could present a total 
filling-in due to microflaring. We have not found filled- 
in by chromospheric emission in the Na I Di and D2 lines 
nor in the Mg I b triplet lines. 

All the activity indicators analysed here indicate that 
the hot component is the active star of the system, con- 
trary to the usual behaviour observed in chromospheri- 
cally active binaries. The variations (in different epochs 
and with the orbital phase) observed in the different in- 
dicators, formed at different heights in the chromosphere, 
are correlated. 

The analysis of several Ti i lines with a positive effect 
in luminosity class indicates that the hot component is 
not a main sequence star but a subgiant or even of higher 
luminosity class, in contrast with the cool component that 
seems to be a dwarf. This result is confirmed by the appli- 
cation of the Wilson-Bappu effect to our Ca 11 K spectra. 
This difference in luminosity class, and therefore in con- 
vcctive zone depth, between the hot and cool components 
could explain why the hot component is the active star of 
the system. 
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